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We compute the dark matter relic densities of neutralinos and axions in a supersymmetric model 
with a gauged anomalous U(l) symmetry. The model is a variant of the USSM (the U(l) extended 
NMSSM), containing an extra U(l) symmetry and an extra singlet in the superpotential respect to the 
MSSM, where gauge invariance is restored by Peccei-Quinn interactions using a Stiickelberg multiplet. 
This approach introduces an axion (Im b) and a saxion (Re b) in the spectrum and generates an axino 
component for the neutralino. The Stiickelberg axion (Im b) develops a physical component (the gauged 
axion) after electroweak symmetry breaking. We classify all the interactions of the Lagrangian and 
perform a complete simulation study of the spectrum, determining the neutralino relic densities using 
micrOMEGAs. We discuss the phenomenological implications of the model analyzing mass values 
for the axion from the milli-eV to the MeV region. These depend sensitively on the value of tan j3. 
The possible scenarios that we analyze are significantly constrained by a combination of WMAP data, 
the exclusion limits from direct axion searches and the veto on late entropy release at the time of 
nucleosynthesis . 
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1 Introduction 



Axions have been studied along the years both as a realistic attempt to solve the strong CP problem 
, to which they are closely related, but also as a possible candidate to answer more 
recent puzzles in cosmology, such as the origin of dark energy, whose presence has found confirmation 
in the study of Type I supernovae 0, 10]. In this second case it has been pointed out that they can 



contribute to the vacuum energy, a possibility that remains realistic if their mass m a - which in this case 
should be ~ 10 -33 eV and smaller - is of electroweak 



111 ] and not of QCD origin. In this case they differ 



significantly from the standard (Peccei-Quinn, PQ) invisible axion. 

According to this scenario, the vacuum misalignment (see 12J, [13J for a discussion in the PQ case) 
induced at the electroweak scale would guarantee that the degree of freedom associated with the axion 
field remains frozen, rolling down very slowly towards the minimum of the non-perturbative instanton 
potential, with m a much smaller than the current Hubble rate. Given the rather tight experimental 
constraints which have significantly affected the parameter space (axion mass and gauge couplings) for 
PQ axions , [nl, the study of these types of fields has also taken into account the possibility to 

evade the current bounds 0, 13- These are summarized both into an upper and a lower bound on the 
size of f a , the axion decay constant, which sets the scale of the misalignment angle 9, defined as the ratio 
of the axion field (a) over the PQ scale vpq {vpq ~ f a ). 

Axion-like particles can be reasonably described by pseudoscalar fields characterized by an enlarged 
parameter space for mass and couplings, with a direct coupling to the gauge fields (of the form aFF) 
whose strength remains unrelated to their mass. They have been at the center of several recent and less 

22, 23], Q, Q, 25]). They are supposed to inherit most of the 
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recent studies (see for instance [19|, [2C ] 
properties of a typical invisible axion - a PQ axion - while acquiring some others which are not allowed 
to it. 

We recall that the axion mass (which in the PQ case is 0(Aq CD / f a ) and the axion coupling to the 
gauge fields are indeed related by the same constant f a . In the PQ case f a (~ 10 10 — 10 12 GeV) makes 
the axion rather light (~ 10~ 3 — 10~ 5 eV) and also very weakly coupled. The same (large) scale plays 
a significant role in establishing the axion as a possible dark matter candidate, contributing significantly 
to the relic densities of cold dark matter. A much smaller value of f a , for instance, would diminish 
significantly the axion contribution to cold dark matter due to the suppression of its abundance (Y x ) 
which depends quadratically on f a . 

It is quite immediate to realize that the gauging of the axionic symmetries by introducing a local 
anomalous U(l) - inherited from an underlying anomalous structure, i.e. a gau ge a nomaly - allows to 



leave the mass and the coupling of the axion to the gauge fields unrelated [26|, l27( , offering a natural 
theoretical justification for the origin of axion-like particles. We just recall that effective low energy 
models incorporating gauged PQ interactions emerge in several string and supergravity constructions, for 
instance in orientifold vacua of string theory and in gauged supergravities (see for instance 28|, |29I]). 

The analysis that we perform in this work has the goal to capture the relevant phenomenological 
features of the axions present in this class of models, extending a previous study presented in a non- 
supersymmetric context [30]. We will be following, as in previous studies, a bottom-up approach. This 
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allows to identify the low energy effective action on the basis of a rather simple operatorial structure 
typical of anomalous abelian models. The theory is then fixed by the condition of gauge invariance of the 
anomalous effective action, amended by operators of dimension-5 (Wess-Zumino or PQ-like terms) which 
appear in the action suppressed by a suitable scale, the Stiickelberg mass (Mst)- 

The introduction of the Stiickelberg multiplet (or axion multiplet), while necessary for the restoration 
of gauge invariance, is in general expected to raise some concerns at cosmological level because of the 
presence, among its components, of a scalar modulus, the saxion. In supersymmetric (ordinary) PQ 
formulations this has a mass of the order of the weak scale or smaller and poses severe problems to the 
standard cosmological scenario. A late time decay of this particle, for instance, could cause an entropy 
release with a low reheating temperature (Trh < 5 MeV) which is unacceptable for nucleosynthesis. Just 
for comparison, we mention that in the case of string moduli, for instance, the interaction of these states 
with the rest of the fields of the low energy spectrum is suppressed by the Planck scale. In turn, this 
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32] in order to enhance the phase 



forces the mass of these states to be quite large (100 TeV or so) 
space for their decay, for a similar reason. 

In our construction the scalar modulus of the axion multiplet acquires a mass of the order of the 
Stiickelberg scale and has sizeable interactions with the other fields of the model, thereby decaying pretty 
fast (~ 10 _23 s). Therefore, smaller values of its mass - in the TeV range - turn out to be compatible with 
the standard scenario for nucleosynthesis. 



1.1 The USSM-A 

Non-supersymmetric versions of the class of models that we are going to analyze have been discussed in 
details in [33j, [26|, |27|] . Recently [34|, [35(, an extension of a specific supersymmetric model, the USSM (the 



U (1) -extended Next-to-Minimal Supersymmetric Standard Model of 36]) has been presented, in which 
the U(l) symmetry is anomalous. This model supports an axion-like particle in its spectrum. It has been 
named the "USSM-A" , to recall both its supersymmetric origin and its anomalous abelian gauge structure. 
It is also close to a similar U(l)' extension of the MSSM (the U(l)' Minimal Supersymmetric Standard 
Model) 37], which supports an axino component among the interaction eigenstates of the neutralino 
sector, but not a gauged axion, due to the structure of the MSSM superpotential. The study of relic 
densities in this model have been performed in 38|, [39fl. In the non-supersymmetric case the identification 
of a physical axion in the spectra of these models has been discussed in detail in (3^], a realization called 
"the Minimal Low Scale Orientifold Model" or MLSOM for short. 

Both in the USSM-A and in the model of 37J], the extra U(l) symmetry takes an anomalous form and 



the violation of gauge invariance requires supersymmetric PQ interactions, with a Stiickelberg supermul- 
tiplet for the restoration of the gauge symmetry. The extra gauge boson of the anomalous £7(1 ) sy mmetry 
is massive and in the Stiickelberg phase, as in previous non-supersymmetric constructions 33l.l26l.l27l]. As 
shown in the case of the MLSOM, axion-like particles appear in the CP-odd spectrum of these theories 
whenever Higgs-axion mixing [33] occurs. For this reason in this work we will be using the term "gauged 
supersymmetric axion" (or axi-Higgs, denoted equivalently as x or Hq) to refer to this state. 

As we have mentioned in the introduction, we will follow a minimal approach. This approach allows 
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to define an effective theory on the basis of 1) an assigned gauge structure (the number of anomalous 
abelian interactions); 2) some conditions of anomaly cancellation and gauge invariance of the effective 
Lagrangian; 3) the choice of a suitable value of the Stiickelberg mass scale characterizing the range in 



which the description of these effective models is compatible with unitarity [40] . As in a previous analysis 



for the LHC in the MLSOM 4l|, we will first stress on the general features of these models, deriving 
the defining conditions for the counterterms which appear in the structure of the effective action, before 
moving to a specific realization with a selected charge assignment. In our simulations we have found that 
the dependence of the results on the choices of the independent charges is, however, extremely mild. In 
this respect the properties that we are able to extrapolate from this class of models - even with a single 
charge assignment - are pretty general and depend quite sensitively only on the choice of the Stiickelberg 
mass Mst and the MSSM Higgs vev ratio tan j3. 

In this section we will focus on the axion/saxion Lagrangian, leaving a general discussion of the various 
contributions to an appendix. It is given by 

^ axion/saxion £-St ~\~ £-WZ (1) 



where Cst is the super symmetric version of the Stiickelberg mass term [42|, while Cy/z denotes the WZ 



counterterms responsible for the axion-like nature of the pseudoscalar b. Specifically 
C st = ^j d A 9{h + bt + V2M st B) 2 



£-wz = — 2 / 



-jg- Tr(gg)b + %L Tr( WW)b 

Mist M S t 



CY -bWtW Y ' a + ^-bW B W B ' a + ^bw Y w B > a 



8{d 2 ) + h.c. , (2) 



M St a M st M st 

where we have denoted with Q the supersymmetric field-strength of SU(3) C , with W the super symmetric 
field-strength of SU(2), with W Y and with W B the supersymmetric field-strength of U(l)y and U(1)b 
respectively. The Lagrangian Cst is invariant under the U(1)b gauge transformations 

6 B B =A + A f 

5 B b = - 2M st k (3) 

where A is an arbitrary chiral superfield. So the scalar component of b, that consists of the saxion and 
the axion field, shifts under a U(1)b gauge transformation. 

The coefficients cj = (cg,cw,cy,cb,cyb) are dimensionless, fixed by the conditions of gauge invariance, 
and are functions of the free charges B{ of the model (as shown below in Eq. ©). Extracting the group 
factors we have 

Abbb Abyy Abby 

CB = "384^ Cy = "l28^ CYB = -Y2^ 

■AbwW A.BGG /.v 
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The coefficients A are denned by the conditions of gauge invariance of the effective action, related to the 
anomalies {U(l)%}, {U (I) B ,U {!%} , {U(l)%, U(1) Y }, {U(1) B , SU(2) 2 } , {U(1) B , SU(3) 2 } . Using the 
conditions of gauge invariance these coefficients assume the form 

Abbb = -3B 3 Hl - 3B 2 Hi (3B L + 18B Q - 7B S ) - 3B Hl (3B 2 + (18B Q - 7B S )B S ) 
+ 3B 3 L + B S (27B 2 Q - 27 B S B Q + 8B 2 S ) 

Abyy = \{-3B L - 9B Q + 75s) 

Abby = 2B Hl (3B L + 9B Q - 5B S ) + (12B Q - 5B S )B S 

Abww = 2^ Bl + 9 ^<3 ~~ 
3 

Abgg = 2 B s- (5) 

We have expressed all the anomaly equations in terms of 4 charges, Bi = (B Bl , B$, Bq, Bl) ordered from 
1 to 4 (left to right), which are defined in Tab. [2j Notice that these charges can be taken as fundamental 
parameters of the model. Their independent variation allows to scan the entire spectra of these models 
with no reference to any specific construction. These relations appear in the anomalous variation (5b) of 
the super symmetric 1-loop effective action of the model, which forces the introduction of supersymmetric 
PQ-like interactions (WZ terms) for its overall vanishing. Formally we have the relation 

SBiB^Suoop + 5 B (ci(Bi))Swz = 0, (6) 

where the anomalous variation can be parameterized by the 4 charges Bi together with the coefficients 
cj(Bi) in front of the WZ counter-terms. In these notations, the uppercase index J runs over all the 5 
mixed-anomaly conditions B s , BY 2 , B 2 Y, BW 2 , BG 2 , ordered from 1 to 5 (left to right). 
Before coming to the definition of the charge assignments we pause for a remark. As we are going to show 
in the next sections, the scalar potential takes a nonlocal form unless all the anomaly coefficients in Eq. [5] 
are zero. Such potential can however be expanded in powers of THeb/Mgt, and as such these contributions 
turn out to be irrelevant if M$t is a very large scale. The situation is rather different if M$t is bound to 
lay around the 1 TeV region, where the potential could actually develope a singularity. In fact, in this 
case, it is in general expected that a singular potential will soon dominate the dynamics of the model. 

We will give the explicit expression of the D-terms for a general choice of the counterterms. The 
function (/) which allows to identify all the charges in terms of the free ones is formally given by 

f(B Q ,B L ,B Hl ,B s ) = (B q ,B Ur ,B Dr ,B l ,B r ,B Hi ,Bh 2 ,B s ). (7) 

These depend only upon the 4 free parameters Bq,Bl, Bh 1 and B$- In our analysis, the charges of 
Eq. have been assigned as 

/(2, 1,-1, 3) = (2,0,-1,1,0,-1,-2,3). (8) 

As we have already mentioned, the dependence of our results on this choice of thi parametric charges 
is very small. Instead, as we will see, the relevant parameters of our analysis turn out to be: 1) the 
anomalous coupling of the gauge boson gs, which controls the decay rate of the saxion and of the axion 
and 2) the Stiickelberg mass. 
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2 Axions, saxions and all orders interactions 



The contributions of the axion and saxion to the total Lagrangian are derived from the combination of 
the Stiickelberg and Wess-Zumino terms. The complete axion/saxion Lagrangian expressed in terms of 
component fields is given by 

£ axion / saxion = ^St £-WZ (9) 

and contains a mixing among the D-terms which is rather peculiar, as we are going to show. The off-shell 
expression of this Lagrangian is given by 



£-axion/ saxion = ^ (<V m & + M st B^) + -d^Reb d^Reb + -Iph^d^h + -^b^d/^b + + ^axion,i 

(10) 

where the expression of C a xion,i is quite lengthy and can be found in the appendix (Eq. ([82]) ). 

The equation of motion for the auxiliary field Fy, can be derived quite immediately and give for the 
F— term of the Stiickelberg field the expression 

Fh ~ ~ 16 Wt G G " 16 Mit WW " 2Wt Y Y " 2Wt * B + 2Wt YB - (U) 
One important feature of the supersymmetric model is that b is a complex field with its real and imaginary 
parts. While Im& may appear in the CP-odd part of the scalar sector and undergoes mixing with the 
Higgs sector, its real part, Reft, the saxion (or scalar axion) before the EW symmetry breaking, has a 
mass exactly equal to the Stiickelberg mass, as expected from supersymmetry. We just recall that in 
the absence of SUSY breaking parameters, the components of the Stiickelberg multiplet form, together 
with the vector multiplet of the anomalous gauge boson, a massive vector multiplet of mass M$t- This 
is composed of the massive anomalous gauge boson, whose mass is given by the Stiickelberg term, the 
massive saxion and a massive Dirac fermion of mass Mst- The fermion is obtained by diagonalizing the 
2-dimensional mass matrix constructed in the basis of Xb - the gaugino from the vector multiplet B - and 
iftb, which is the axino of the Stiickelberg multiplet. The diagonalization of this matrix trivially gives two 
Weyl eigenstates of the same mass Mst, which can be assembled into a single massive Dirac fermion of 
the same mass. Notice that in this re- identification of the degrees of freedom contained in b and in the 
vector multiplet B, Imb takes the role of a Nambu-Goldstone mode and can be gauged away. 

The saxion has typical interactions of the form RebFiFj, with the gauge fields which have mixed- 
anomalies with U(1)b, beside non-polynomial interactions with the remaining fields of the theory. As we 
are going to elaborate, this features shows up because of the presence of terms consisting of the product 
of two D fields and the saxion. To clarify this point, we recall that the general Lagrangian contains a 
supersymmetric Wess-Zumino term of the form 

1 / iAn C Y ,-,„y w y a 



Cwz,y = -~ / (TO -±-hW Y a W Y a 8(6') + h.c. (12) 
2 J M S t 



which gives, after the expansion in components, a term proportional to 



^-RebD Y D Y . (13) 

Mst 
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This kind of terms, once that the equations of motion (EOM) of the auxiliary fields D are calculated and 
substituted back into the lagrangian, give the non-polynomial form of the potential. Furthermore, from 
the WZ term corresponding to the anomaly BBY (which is the term proportional to cyb), we get a term 
proportional to Reb DyD b so that the EOM for the abelian D fields are coupled. 

The derivation of such equations involves all the terms of the Lagrangian discussed in appendix [XJ We 
obtain 



Db,os 



1 



12 + 12y/2Reb(c B + c Y )/M i 



st 



6Reb 2 (c 2 YB -Ac Y c B )/M 2 t 



} _CB_ 

'Mat 



V2 + 2 



CY 

M Si 



Reft 



M 2 st 



Reb 



(-3i\ B ^b + h.c.) - 12M st Reb - 12v / 2c y Re6 2 + 



12g B M + v^^g-ReM (B S S^S + B Hi H\Hx + B H2 H\h 2 + B Dr D r D r + B Ur U r U r + 



B Q Q ] Q + B R R^R + B L Vi) + 3\/2^-(UyVb + h.c)+ 

Mst 

V2j^-g Y Reb( - 3H\h x + 3B H2 H ] 2 H 2 + 2D R D R - 4U R U R + Q+Q + 6i? f R - 3DL) } 



Dvos 



12 + 12y/2Reb(c B + c Y )/M st - 6Reb 2 (c YB - ic Y c B )/M 2 t 

cyb 



3 -TPT Re6(iAyVb + h.c.) + 3\/2 -^-(iXyipb + h.c.) - 6\/2 cy B Re6 2 + 
M st M S t 

CYB 



{ 

6V2 Reft g B {Bs&S + B Hi H\H x + B H2 H\h 2 + B Dr D r D r + B Ur U r U r + 

Bq&Q + B R R)R + B L V~L) - 6^f- (V2 + 2-^- Reb B ) {i\ Y ^ h + h.c.)+ 

■MSt V M St J 

2g B (l + v^^-Reft^ {-SHIR! + 3B H2 H ] 2 H 2 + 2D R D R - 4U R U R + Q ] Q + 6^ - 3tf~L)}, 

(14) 

showing that their on-shell expressions are characterized by the appearance of the saxion field in a non- 
polynomial form. The presence of the Stiickelberg mass allows to perform an expansion of these terms to 
all orders in Reb /Mst- The first terms of the series expansion are given by 

1 



12 + \2^/2Reb{c B + cy)/M st - QReb 2 {c 2 YB - Ac Y c B )/M 2 t 

1 Re6 Re6 1 , Refe 2 1 Re6 2 

Cy — Cr — 

12 6^/2M st 6V2M, 



st 



1 2 Reb 2 
+ 6 Cy M 2 " 



st 



+ 24# B 



1 2 Reb 2 1 _ Reb ) T , 1f2 

Wt + t CB W t ' 



H CnCy „ 

6 M 2 t 



+ 0(Reb 2 /M' St ) 
(15) 



We present a list of the vertices to leading order in 1/M$t in Fig. [TJ Additional vertices (not shown) 
come with n insertions of Re b and a suppression by higher powers (2n) of Mst- Some considerations 
are in order concerning the allowed values of Mst- A very large Stiickelberg mass, in principle, would 
be sufficient to guarantee that the effect of reheating - caused by the decay of the saxion - takes place 
well above the temperature of nucleosynthesis (see for instance the discussion in 32J) thereby avoiding 
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Re6 




Rc6 




Reft 



Re6 




Reft ■ 



Rci> 



Qt 



Rci) 




Rci 



Figure 1: Saxion interactions to lowest order in 1/Mst- An infinite number of additional higher order interactions 
(in powers of 1/Mgt) are generated by the insertion on these vertices of n powers of saxion lines. We use the double 
line notation for Majorana particles. 



the problem of a possible late entropy release at that time. In this case one can essentially neglect the 
saxion from the low energy spectrum. For moduli of string origin the required mass value (~ 100 TeV), 
much larger than in our case, is justified by the suppressed gravitational interaction of the modulus with 
the rest of the matter fields and works as an enhancing factor for its decay. In our case, instead, such a 
suppression is absent and a fast decay of the saxion is guaranteed already by a Stiickelberg mass around 
the TeV scale. 
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3 The scalar potential and the saxion 



As we have mentioned, in this model there are three scalar fields which take a vev, H±,H 2 and S, the 
scalar components of the scalar superfield S. We assume that also the saxion field gets a vev, v b . The 
scalar potential is composed by contributions coming from the D-terms, F-terms and scalar mass terms. 
Expanding up to 0(Reb/Mst) we find 



V 



|AFi-iT 2 | 2 + |AS| 2 (|i7i| 2 + |tf2| 2 ) 



m 2 |#i| 2 + m 2 2 \H 2 \ 2 + m 2 s \S\ 2 + (a A Sfli • H 2 + h.c.)- 



±g Y (fffffi - H\H 2 ) 2 - X -g 2 B (B Hl R\R x + B H2 H\h 2 + B S &S)' 



l -g 2 Y (HlH! - 


1 2 

8 92 


(hIh, - 


Reft 3 


~c B M S t 
[ V2 


Re6 2 


~ m Reb , 

2 


Re&j 


g B M st ( 



+ 



+ 9Y 



(cy + C B ) C Y B 



+ 



2M st V 1 1 2 V_ 
c B 5B (b Hi hIH! + B H2 H\h 2 + S^s) + ^ , T 



+ 



C-W92 



H\Hi - h\h 2 I + 



CBS! 



v/2M, 



> 32 v / 2M 5t 
Bh 1 h\Hi + B H2 H\h 2 + BgSte) 2 + 



(HlH! - HlH 2 y + A\H\H 2 \ 



+ 



cy9 y 
4V2M St 

CYB 9Y9B 

2V2M St 

To proceed with the analysis of this potential we introduce the following parameterizations 



(b Hi hIh! + b H2 h\h 2 + b s s^s) [h\H\ - h\h 2 ) j. 



(16) 



1 / Rei7° + i ImHf 
V2 \ReH{ + i Imfff 



1 / Re# 2 + + i Im#+\ 
72 \ Retf £ + i ImflJ J 



1 



-^(ReS + i ImS). 
v2 



expanded around the vevs of the Higgs fields and of the saxion as 



(Hi) 





tan/? = — w) = — f (Reo) = 

V2 



• (17) 



The scalar mass parameters can be expressed in terms of the remaining parameters of the theory using 
the minimization conditions for the scalar potential. In particular, taking a derivative of the potential 
with respect to the saxion field we get the relation 



dV 



<mvh -v b rnl eb - v b M 2 t + -g B M st B Hl v 2 + -g B M st B H2 v 2 + -g B MstB s v 2 s (18) 

where we have neglected all the contributions suppressed by the Stiickelberg mass. We can use this relation 
as a necessary condition in order to express ni 2 ^ eb in terms of the vevs and of the other parameters of 
the scalar potential. A numerical analysis of the Hessian at this point, for the selected parameters of the 
model used in our simulations, shows that indeed this extremal point indeed corresponds to a minimum. 

We will try to highlight the most interesting features of these types of models and the implications 
for the axion, which are all connected to the properties of the vacuum of these theories below the scale of 
SUSY breaking and at the scales of the electroweak and QCD phase transitions. 
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4 Saxion decay modes 



Having summarized the basic structure of the model, we now turn to describe the leading contributions 
to the 2-body decays of the saxion . The goal of this analysis is to ensure that the decay rate of the saxion 
is such that it occurs fast enough in order not to interfere with the nucleosynthesys. 

We will compute its decay rate by considering the worst possible scenario, i.e. by assuming that this 
decay occurs around the SUSY breaking scale, or temperature T around 1 TeV. At this temperature, the 
decays of the saxion are parameterized by the typical SUSY breaking scales such as Mj, and My, both 
of 0{M susy ). The model is in a symmetric electroweak phase {M susy > v), which justifies the use of the 
interaction eigenstates (rather than the mass eigenstates) for the description of the final decay products. 

The relevant interactions for the saxion decay are described by the general Lagrangian 

g B M S t B Hi H{h± + g B M S t B H2 H\h 2 + g B M st B S S^S+ 

3 3 3 

g B M S t B Q Q)Qj + gsMst B Dr DrjDrj + 9 B M st B Ur ^ U^Urj+ 
j=i j=i j=i 

g B M st B L Y l\Li + g B Mst B R ^ ^Rl 



sn.T.irtn.dp.r.. ^ 



-l- 



l=e,fi,r l=e,[j,,T 
C B I, i t \ T \ ■ C YB 



(19) 



: — /| (Ab ipb + ipb) - i—f= (ip b Ay + iph Ay) 

They involve CP-even and CP-odd massless scalars, the extra singlet scalar S, the squarks, the sleptons 
and the gauginos ipb, Ay. We compute the total decay rate into fermions, squarks and sleptons, and Higgs 
scalars. The left-handed doublets of the squarks and the sleptons are defined as Qj and Li respectively, 
while the right-handed singlets are Urj, -Drj and R[, with j,l labeling the fermion families. 

• Decays into fermions 

Assuming that Mb ~ My are slightly less than 1 TeV, the decay rates of the saxion into one gaugino and 
one axino are 

r(Re^W b H4^(l-4-f-f ! . 

r(Re^A,* b )=4 B ^(l-4-| : ) S/2 (20) 

with the expressions of the coefficients c B and cy B determining the couplings given explicitly in Eq. @. 
Notice that these rates are large due to the linear dependence on MR e {,(= M$t)- 

• Decays into Squarks and Sleptons 

In this channel we consider, for simplicity, the decay only into squarks and sleptons of the same type. 
Even in this case we are assuming that the masses of the squarks and of the sleptons are all equal and 
slightly below 1 TeV. The decay rate into the i-type sfermion is given by 

r ( Re6 ^) = ifSt (21) 
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where the kinematic function A is, in general, defined as A = (m? + M| - M| e J - AMfMf (h ere with 
Mj = Mj), and the couplings gi, in the various cases, are defined as 

N c cjj r R-handed singlet u-type squark, 

N c cd r R-handed singlet quark d-type squark, 

N c cq l L-handed doublet squark, (22) 

cr R-handed singlet slepton e, /2, f, 

cl L-handed doublet slepton. 

Here N c = 3 is the color factor and the various couplings are given as 

cd r = 7^9b Mst B Dr , c Ur = ^g B M S t B Ur , c Ql = -^g B M st B Q , 
c L = -\qb M st B L , c R = X -g B M st B R . (23) 
• Decays into massless scalars 

The decay rate into particles of the Higgs sector that we denote generically with hi = Rei^i , Imiii , . . . 
is given by 



a? I Ml 



1/2 



TCReb^hihi) = l - 1 - 4^f^ , (24) 

V l) 32vrM Rcfc I Ml eb K ' 



where the couplings Si are defined as 



Si 



ch 1 Hi Higgs doublet 

ch 2 H2 Higgs doublet (25) 



1-2 

cs S Higgs singlet 



and the coefficients ch 1 , ch 2 > cs are 



ch x =-j9B M st B Hl , c H2 = -j9B M st B H2 , c s = ~^9b M st B s . (26) 
The total decay rate is obtained by summing over all the decay modes 

r tot = r (Re b -> A^ b ) + r (Re b -> Ay ^ b ) + r (Re b fji^j + ^ r (Re b -> . (27) 

i j 

All the decay rates depend upon the value of the extra Ub{X) coupling g B , the Stiickelberg mass Mst and 
the SUSY breaking scale M susy . 

The total decay rate and the lifetime of the saxion are shown in Fig. (|2|), with the saxion mass Mn e (, 
given by the Stiickelberg scale (M^ = Mst) around 1.4 TeV, and with all the squarks and the sleptons 
in the final state taken of a mass of 700 GeV. All the particles of the Higgs sector are considered to 
be massless. For g B = 0.01 we obtain a saxion whose decay rate is around 60 MeV if its mass is 1.7 
TeV, and which decays rather quickly, since its lifetime is about 10 -23 seconds. The lifetime decreases 
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Figure 2: Total decay rate and lifetime of the saxion for different values of gB as a function of the Stuckelberg mass. 



quite significantly as we increase the gauge coupling of the anomalous gauge symmetry. For instance, for 
gs = 0.1 it decreases to ~ 10 -24 sec, since the phase space for the decay is considerably enhanced. 

We conclude that the saxion decays sufficiently fast and does not generate any late entropy release at 
the time of nucleosynthesis. Obviously, this scenario remains valid for all values of the Stuckelberg mass 
above the 1 TeV value. Therefore, in the analysis of the evolution of the contributions to the total energy 
density (p) of the universe, either due to matter (p m ) or to radiation (pr), at temperatures T < 2 TeV, 
the contribution coming from the saxion is entirely accounted for by pr. 

At this point, having cleared the way of any possible obstruction due to the presence of moduli at the 
low energy stage (T < 2 TeV) of the evolution of our model, we are ready to discuss the relevant features 
of the Stuckelberg field. In particular, we will discuss the appearance of a physical axion, the physical 
component of the Stuckelberg, at the electroweak scale. This is extracted from the CP-odd sector and 
generated by the mechanism of vacuum misalignment taking place at the same scale. In particular, the 
discussion serves to illustrate how a flat - but physical - direction might be singled out from the vacuum 
manifold, acquiring a small curvature at the electroweak phase transition. 



5 The flat direction of the physical axion from Higgs-axion mixing 

In 3J] , 35f] we have presented in some detail an approximate procedure in order to identify in the CP-odd 
sector one state that inherits axion-like interactions. The approach did not require the explicit expressions 
of the curvature terms in the CP-odd part of the supersymmetric potential, which are instead needed 
in the discussion of the angle of misalignment. Here we are going to extend this analysis by giving the 
explicit parameterization of these additional terms. The determination of the angle of misalignment and 
its parameterization in terms of the physical axion is based on an extension of the method presented in 



301 ] . We are going to illustrate this point starting, for simplicity, from the non-super symmetric case and 



then moving to the supersymmetric one. 
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5.1 The non-supersymmetric case 

In the non-supersymmetric case the scalar sector contains two Higgs doublets Vpq(H u , H d ) plus one extra 
contribution (a PQ breaking potential), denoted as V^^(H U , H d ,b), which mixes the Higgs sector with 
the Stiickelberg axion 6, 

V = V PQ (H U , H d ) + V^(H U , H d , b). (28) 

The mixing induced in the CP-odd sector determines the presence of a linear combination of the Stiick- 
elberg field b and of the Goldstones of the CP-odd sector, called x-> which is characterized by an almost 
flat direction, whose curvature is controlled by the strength of the extra potential Vpq(H u , H d ) is 

the ordinary potential of 2 Higgs doublets, 

V PQ = nlKlK u + &H\H A + \ uu {HlH u ) 2 + X dd (H d H d ) 2 - 2X ud (HlH u )(H d H d ) + 2X' ud \H^r 2 H d \ 2 

(29) 

Concerning the VrfH contribution to the total potential, its structure is inferred just on the basis of gauge 
invariance and given by 

V.. =\ Q {HlH ie ~ i9B(BH ^ BH ^^) + \ x {HlH x e~ i9B{BH *- BH ^^) 2 + 

X 2 (HlH 2 )(HlH ie ~ igB{BH ^ BH ^^) + A 3 (F 1 t ^ 1 )(if 2 t J ff 1 e" i9B(i3H ^ BHl) ^) + h.c. (30) 

These terms are the only ones allowed by the symmetry of the model and are parameterized by one 
dimensionful (Ao = Xqv) and three dimensionless constants (Ai, A2,Aa). 

The CP-odd sector is then spanned by the three fields (Imiii , Imil2 j b) , with the potential Vpq a 
function only of H\ and H 2 . After electroweak symmetry breaking, due to Higgs-axion mixing, b can be 
written as a linear combination of a physical axion and of an extra component. The latter is a linear 
combination of the two Goldstone modes of the total potential (Vpq + V^), denoted as Gq,Gq. The 
physical axion, x, i.e. the component of b which is not proportional to the two Goldstones, can be identified 
using the rotation matrix O x which relates interaction and mass eigenstates in the CP-odd sector 



(31) 



which takes the form 

b = O^Gl + O^G 2 + 0^ X . (32) 

X inherits WZ interactions from b via Eq. (|32[) . once this is introduced into the WZ counterterms. 

From an explicit computation one finds that Of 3 = 0, 0£ 3 ~ O(l) and 0% 3 ~ v/Mst- The Goldstones 
of the two neutral gauge bosons, Gz,Gz' are linear combinations of Gq and Gq and can be extracted 
from the bilinear mixings after an expansion around the broken electroweak vacuum. Then, the entire 
CP-odd sector can be spanned by the basis (Gz,Gz',x)- The presence of an extra degree of freedom in 
this sector has been established in 
point for clarity. 









k 


= o x 








\ b J 



331 ] using a simple counting of the degrees of freedom. We review this 
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There are 9 degrees of freedom in the set (Ay, W%, B, ImiTi, Imi^), where B is the massive Stiickelberg 
gauge vector field, before electroweak symmetry breaking, as well as 9 in the set (Ay, Z, Z', x), which is 
generated after the breaking. The direction determining the gauged axion \ is then physical but flat, 
in the absence of an extra potential which may depend explicitly on b. The potential V' is responsible 
for giving a small mass for \ an d can be used to parameterize the mechanism of vacuum misalignment 
originated at the electroweak scale. 

One can explore the structure of this potential and, in particular, investigate its periodicity. The 
phase of the potential is indeed parameterized by the ratio x/ a x 13Q(] 



V if4 = 4v i y 2 {^vf + X 3 vj + A ) cos (j^j + 2Ai«?u| cos (zf^J 



(33) 



with a mass for the physical axion \ given by 



in 



2 _ 2v 1 v 2 



a 2 

x 



(X v 2 + X 2 v\ + \ z v 2 2 + 4A 1 wiw 2 ) ps Xefjv 2 , (34) 



with a x ~ 0(y). The size of this expression is the result of two factors which appear in Eq. (]34p : the 
size of the potential, parameterized by (Ao, Ai, A2, A3), and the electroweak vevs of the two Higgses. The 
appearance of \ i n Eq. (|33p - in the phase of the extra potential - shows explicitly that the angle of 
misalignment is entirely described by this field. The angle is defined as 

0(x) = (35) 



where 



a = 2v 1 v 2 Mst , 3g . 

lgl(B H2 - B Hx Yv\vl + 2Ml t (vj + v 2 ) 



is the new dimensionful constant which takes the same role of the scale f a of the PQ case (9(x) = a/f a ). 
The potential is characterized by a small strength ~ A e j^w 4 , and for this reason one can think of x as a 
pseudo Nambu-Goldstone mode of the theory. 

At this stage, it is important to realize that the size of the extra potential is significant in order 
to establish whether the degree of freedom associated to the axion field remains frozen or not at the 
electroweak scale. For instance, if X e jf is associated to electroweak instantons (X e ff ~ Aj ns t), then m x 
is very suppressed (see the discussion in Sec. I5.3f) and far smaller than the corresponding Hubble rate at 
the electroweak scale 

which is about 10 -5 eV. In the expression above g* Ti is the number of effective massless degrees of freedom 
of the model at a given temperature (T), while Mp denotes the Planck mass. We recall that the condition 

m x (T) ~ 3H(T) (38) 

which ensures the presence of oscillations and determines implicitly the oscillation temperature Tj, is 
indeed impossible to satisfy if the misalignment that generates the value of m x at the electroweak scale 
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is assumed of being of instanton origin (see the discussion in Appendix [Cj . This implies that the degree 
of freedom associated to this physical axion would be essentially frozen at the electroweak scale, and 
the oscillations could take place at a later stage in the early universe, only around the QCD hadron 
transition. Instead, a more sizeable potential, providing an axion mass larger than 10 -5 eV, would allow 
such oscillations. For an axion mass around 1 MeV oscillations indeed occur, but are damped by the 
particle decay, given that its lifetime (77 ~ 10~ 4 sec) is much larger than period of their oscillation ( 
Tosc ~ 1CP 13 sec). This discussion is going to be expanded to the supersymmetric case. 

5.2 Supersymmetry and the angle of misalignment 

In the supersymmetric case the situation is analogous, in the sense that the physical direction x can be 
identified by the same criteria. The superpotential that we are considering allows the presence of one extra 
degree of freedom, given by ImS", to appear in the CP-odd sector besides the states (Im Hi , Im H2 , Im b) , 
already present in the non-supersymmetric case. 

From the supersymmetric potential V in (|16p . we identify two massless states, that we call G\ and 
G2, and a massive eigenstate, called H®. G\ and G° do not coincide with the true Goldstones of the 
model, as in the previous case, since the V potential does not include any contribution involving Im&. 
The correct neutral Goldstone modes are extracted from the derivative couplings between the CP-odd 
scalar fields and the neutral gauge bosons present in the Lagrangian. The physical axion is then identified 
as the massless direction which is orthogonal to the subspace spanned by (Gz,Gz' , H®). This state is 
called Hq = x and is given by the linear combination 



It is important to remark that this state is not constructed, at least at this stage, from the matrix O x , 
since the projection of Im& on x would be zero, if the matrix O x were derived just from the potential V 
since it does not depend on Im b. 

Also in this case, the identification of the Goldstones of the two neutral massive gauge bosons {Z and Z') 
is obtained by looking at the bilinear mixing terms; these appear in the Lagrangian once this is rewritten 
in the physical basis (in the form MzZdGz, and Mz'Z'dGz 1 )- Then one can immediately figure out 
that the linear basis spanning the entire CP-odd sector can be completed by the addition of an extra, 
orthogonal state % {Gz,Gz',H2,x)- The new entry parameterizes a massless but physical direction in 
this sector. Once Im& is re-expressed in terms of the physical axion x an d of the Goldstone modes Gz, Gz i 
of the massive gauge bosons, x wm inherit from Im6 axion-like interactions and will be promoted to a 
generalized PQ axion. 

At this point, having identified this flat but physical direction of the potential in the CP-odd sector, 
one can ask the obvious question whether the same potential can acquire a curvature. These effects are 
indeed parameterized by the strength (\ e ff) of the potential V (the "extra potential") that we are going 
to identify below, and which remains a free parameter in the theory. 



X = [M st vivl ImflJ + M st vjv 2 Imtf 2 ° - M st v 2 v s ImS - B s g B (v 2 v 2 s + v\v\ 



< 2 2 )Im&] 



x 




(39) 
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One special comment is deserved by vs, the vev of the scalar singlet, which is new compared to 
the standard MSSM scenario and which is part of the scalar potential. We recall that this new scale 
is essentially bound by the condition Xvs ~ \i ~ 10 2 - 10 3 GeV (see Eq. [78]). This defines the typical 
range for the \i term, which sets the scale of the interaction for the two Higgs doublets in supersymmetric 
theories. 

In our case we are allowed to parameterize this new non-perturbative contribution (V) to the potential, 
as discussed in the previous section, in a rather straightforward way, by classifying all the phase-dependent 
operators which can be constructed using the fundamental fields of the model. In analogy to the non- 



supersymmetric case (the MLSOM) 33] we rely only on gauge invariance as a guiding principle to identify 
them. These include, in particular, a dependence of V' , again in the form of a phase factor, from the 
Stiickelberg field Im b. 

The contributions appearing in V don't need to be given necessarily in a supersymmetric form, since 
we are assuming that supersymmetry is already broken at the scale at which they appear (v < M susy ). 
They are parameterized in the form 

V = V) Vi (40) 



=1 



where 



-iig B B s - r linb 



V^a^e +h.c. 



Imb 



V 2 = e l9BBs iM St ( a2Hl . H2S 2 + hHlHtS + b 3 H%H 2 S + b 4 S^ S 2 + dtS) + h.c. 



V 3 = e i9B2Bs ^ b St (asHlHxS 2 + a^H^S 2 + a 5 S^S 3 + Cl S 2 ) + h.c. 



, 2 ig B 2B s ^- 



V 4 = a 6 (ifi • H 2 ) z e yB ^^ + h.c. 
V 5 = hS 3 e- i9B3Bs ^k+h.c. 

V 6 = C2 Ht ■ H 2 e 9B s 2M st + h.c. (41) 



In the expressions above we have grouped together terms that share the same phase factor. Notice that the 
parameters a^, bj, cj~ and d\ carry different mass dimensions. For these reasons they can be parameterized 
by suitable powers of the SUSY breaking mass M susy times \ e f f . We explicitly obtain the estimates 

a i ~ ^eff bj ~ X e ff M susy C k ~ X e ff M 2 usy di ~ X e ff Ms Usy . ( 42 ) 



If we introduce any of the terms in Eq. (|4ip . and recompute the CP-odd mass matrix using the new 
potential (V + V), this gets modified, but we still find two massless eigenstates corresponding to the 
neutral Goldstone modes, which also in this case we call Gq and Gq. They can be expressed as linear 
combinations of the neutral Goldstone states coming from the derivative couplings between the gauge 
bosons and the CP-odd Higgs fields. An important point to remark is that these states (Goldstone 
modes) do not depend on the parameters of the Peccei-Quinn breaking potential, as we expect, since the 
presence of this extra potential doesn't affect the bilinear derivative couplings through which they are 
identified. 
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In the basis (ImHl , ImH^ , ImS, Imb) they are given by 
^ = (^,0,0 



v 



2 2 

Gq = = <j gBB s ^,9BB s ^,-g B B s v s ,Mst } ■ (43) 



s u s> 



5.3 The strength of the potential and A e // 

One important comment concerns the possible size of the axion mass m x induced by V at the electroweak 
scale. In this respect we will take into account two basic possibilities. A first possibility that we will explore 
is to assume that the axion mass is PQ-like, in the milli-eV region; as a second possibility we will select an 
axion mass around the MeV region. These choices cover a region of parameter space that has never been 
analyzed in these types of models, while a study of the GeV region for the axion mass has been addressed 
before in 4l|. These choices have to be confronted with constraints coming from a) direct axion searches, 
b) nucleosynthesys constraints and c) constraints on the relic densities from WMAP data. 

A PQ-like axion is bound to emerge in the spectrum of the theory if the potential V is strongly 
suppressed and the real mechanism of misalignment which determines its mass is the one taking place 
at the QCD transition. The value of A e /y, under these assumptions, should be truly small and one 
way to achieve this would be to attribute its origin to electroweak instantons. Using the numerical 
relations for the electromagnetic (a) and weak couplings (aw), l/a(Mz) = 128 and a\y = a/sm 2 9\y 
with sin 2 6\y(Mz) = 0.23 on the Z mass (aw{Mz) = 0.034), the exponential suppression of the extra 
potential is controlled by X e ff ~ e~ 185 = \i ns t = 4.5 x 10~ 81 . This corresponds to a mass for the axion 
given by m x ~ \/^eff v ~ 10~ 29 eV. This mass would be obviously redefined at the QCD epoch. 

As we have briefly mentioned in the introduction, mass values of the axion field around 10 -33 eV (for 
global U(iys or of PQ type) and with a spontaneous breaking scale f a ~ 10 18 eV have been considered 
as a possible origin of a cosmological constant A 4 ~ (lO _3 eV) [ll|. In such models the misalignment is 
purely of electroweak origin and connected to electroweak instantons. Oscillations of fields of such a mass 
would not take place even at the current cosmological time. 

Instead, for an axion of a mass in the MeV region, the value of \ e ff is larger (~ 10~ 12 ) and will 
be estimated below. In this case the effect of vacuum misalignment at the QCD scale is irrelevant in 
determining the mass of this particle. A more massive axion, in fact, decays at a much faster rate than 
a very light one and the usual picture typical of a long-lived PQ-like axion, in this specific case, simply 
does not apply. 

In order to characterize in more detail the potential in Eq. (|40p . we proceed with a careful analysis 
of the field dependence of the phase factors in the exponentials, that we expect to be written exclusively 
in terms of the physical fields of the CP odd sector, Hq, and the axion x (x = ^o)- I n ^ct, this is the 
analogous (and a generalization) of what found in the previous section (see Eq. ([33]) ). where the periodicity 
has been shown to depend only on the axion \. For this purpose we use the following parameterization 
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of the fields 



H\{x) = ± ( p \ {x) + Vl ) e »iM H f(x) = _^pf ( x )e^W 



S(x) = ^(p 5 (x) + ^)e^) 



x/2 



(44) 



and select just some of the in Eq. (j4Tj) in order to illustrate the general behaviour. 

For instance, if we consider only the V\ term we get the corresponding symmetric mass matrix for the 
total potential V + V\ , with V defined in Eq. (fl~6j) , 



/ V2VS 
1 VI 



M, 



odd 



V2 



V2 



V-2 







a-\ Mst 



\ 



_ „2 r2.,4 
n /nOl 9B D S V S 



(45) 



expressed in the basis (<£;[, <3?2' ^s, Imb). From this matrix we get two null eigenvalues corresponding to 
the neutral Goldstones and two eigenvalues which correspond to the masses of the two CP-odd states 
and Hq. In this specific case they take the form 

1 



A ± VA 2 - B\ 



2M S tv 1 v 2 vs 

A = Aa lVl v 2 vl (4M| f + g%B 2 s v 2 s ) + V2a x M 2 st (vfv% + v 2 v 2 s ) 



B = lQV2a 1 a\Ml t v 1 v 2 v% (±v 2 M 2 st + g 2 B B\ 



2„,2 



[V, V. 



l u 2 



2 2 

V Vq 



(46) 



In the limit of a vanishing a\ (~ ^eff) we obtain a massless state corresponding to Hq (x) an d a massive 
one corresponding to Hq. In fact, expanding the expressions above up to first order in a±, which is a very 
small parameter due to (|42l) . we obtain for the two eigenvalues the approximate forms 



H 4 



v 2 v 2 v 2 



>'2 



v 2 v% + vtv 



J s 



2„ ,2 • 



l u 2 



m H5 



Aax4 [^ 2 M 2 t + g 2 B B 2 (v 2 v 2 +v 2 v 2 )] 



M 



st 



2„,2 



v 2 vl) 



(47) 



These relations show that indeed m H 5 is 0(X e ffv) while m H A is 0(v). 

Moving to the analysis of the phase factor of the same term (Vi), the linear combination of fields that 
appears in the exponential factor is given by the expression 



n = 



4$ 5 (x) 2g B B s l™b{x) 



vs 



M st 



(48) 
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Figure 3: Illustration of the two misalignments at the electroweak (upper figure) and at the QCD phase transitions 
(lower figure) for a PQ-like axion (not to scale). 



We rotate this linear combination on the physical basis (Gz,Gz, Hq, Hq) using the rotation matrix O x . 
After the rotation we can re-express the angle of misalignment as a linear combination of the physical 
states of the CP-odd sector in the form 

a_H% + Hl (49) 



This linear combination will appear in all the operatorial terms included in V and is a generalization of 
Eq. ()35|) . with a H 4 and a H s defining, separately, the scales of the two angular contributions to the total 
phase. It is not difficult to show that the periodicity of the potential depends predominantly on Hq. This 
can be easily seen by analyzing the size of a H 4 and cr H 5. In fact, expanding to first order in a\ we get 

1 _ A Vl v 2 8V2vy 2 vl [4v 2 M 2 st + g|g| (v 2 v 2 s + v 2 v 2 )] a, 

2 ; „.2„,2\5/2 ~ + ° ^ 



a H$ V2 sgn B s ^v 2 v 2 + v 2 v\ a x M 2 t (v*v 2 + v 2 v 2 ) sgn (B s 

1 8V2v 2 v 2 v 4 s [Av 2 M 2 t + g 2 B B 2 (v 2 v 2 + v 2 v 2 )] ai 



<r H $> a x M 2 t (v 2 v 2 + v 2 v 2 f' 2 sgn B s 



+ O (af) . (50) 



with a\ being proportional to the SUSY breaking scale M susy . A more careful look at the structure of 
these two scales shows that a H 4 ~ vs (vs = 400 GeV in our case) while a H 5 ~ M susy /X e ff. Clearly, 
a H 5 cf H 4, but the dependence of the extra potential on \ is clearly affected by the different possible 
sizes of X e ff- For an instanton generated potential (X e ff ~ \nst) the direction of x is essentially flat 
and a H 5 turns out to be very large. In turn, this implies that the dependence of the potential V\ on Hq, 
which takes place exclusively through the exponential, is negligible, being essentially controlled by Hq 
H%/v) with 



V^X^fV 4 COS(0!). (51) 

We may conclude, indeed, that in this case the effect of misalignment on x, generated at the electroweak 
scale, can be neglected. This feature is shown on the left panel of Fig. HI where we plot V'(Hq,x)- It is 
immediately clear from these plots that for \ e ff ~ ^inst the only periodicity of the extra potential is in 
the variable Hq (left panel), due to the flatness of the Hq direction. For a more sizeable potential, with 
X e fj ~ 10~ 12 , the curvature generated in \ is responsible for giving a mass to the axion in the MeV range 
(Fig. HI right panel). This result is generic for all the terms. One can draw some conclusions regarding 
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Figure 4: Shape of the extra potential V at the electroweak scale in the CP-odd sector in the (Hq,x = Hq) plane. 
X is an almost flat direction for a strength induced by the instanton vacuum at the electroweak scale (left panel) 
and acquires a curvature for an axion mass in the MeV region (curvature in the \ direction, right panel). 

the role played by the exponential phase and compare the supersymmetric with the non-supersymmetric 
case. In the non-supersymmetric case the periodicity of the potential is controlled by the weak scale (v), 
and is expressed directly in terms of the physical component of b (which is a real field). The size of the 
potential, in this case, is of order \ e ffV 



and therefore very small, while the periodicity shows that the amplitude of the axion field is x ~ 0(v). 
In the supersymmetric case, more generally, we obtain for a generic component Vi 



from which it is clear that the curvature in the axion field is controlled by the parameter X e ff- I n the 
supersymmetric case we can think of the periodicity in Eq. (|53p as essentially controlled by the massive CP- 
odd Higgs Hq, with a period which is 0(irvs), with superimposed a second periodicity of 0(7rM susy /X e ff) 
(with M susy /X e ff 3> vs) in the perpendicular direction (x)- We conclude that the actual structure of 
the complete (V + V) potential indeed guarantees the presence in the spectrum of a physical and light 
pseudoscalar field. This analysis holds, in principle, for an axion of any mass, although we do not explicitly 
study an axion whose mass goes beyond the MeV region. To have an axion which is long-lived, the true 
discriminant of our study is the axion mass, and for this reason we are going to present a study of the 
decay rates of this particle keeping the mass as a free parameter varying in the milli-eV - MeV interval. 

6 Decay of a gauged supersymmetric axion 

In this section we compute the decay rate of the axion of the supersymmetric model into two-photons, 
mediated both by the direct PQ interaction and by the fermion loop, which are shown in Fig. [5l keeping 
the axion mass as a free parameter. Denoting with N c (f) the color factor for a fermion specie, and 




(52) 




(53) 
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Figure 5: Contributions to the axi-Higgs decay \ — > 77. 




introducing the function 1777(7/), a function of the mass of the fermions circulating in the loop with 

1 



t = Amj/m 2 , r/(r) = arctan 2 — -== pf x = J 1 — ( J , (54) 



the WZ interaction in Fig. [5] is given by 



where g^ is the coupling, defined via the relations 



M wz(x -> 77) = 4g*e[iJ,, v , ki,ki], (55) 



x _ 9bB s {Ac Y g\ + c w g Y ) / v 2 v 2 s + v\v\ , , 

9r n~ 16g 2 M st V 4M 2 t v 2 + g 2 B B 2 (v 2 v 2 + v 2 v 2 ) 



obtained from the rotation of the WZ vertices on the physical basis (we will comment in more detail on 
the size of this coupling in the next section). The massless contribution to the decay rate coming from 
the WZ counterterm xF^F^ is given by 

Hi 

IWX^77) = ^(9* 7 ) 2 . (57) 

Combining also in this case the tree level decay with the 1-loop amplitude, we obtain for x ~ * 77 the 
amplitude 

M^{ X ^ll) = M^ z + Mf. (58) 
The second amplitude in Fig. is mediated by the triangle loops and is given by the expression 

Mf{x -> 77) = ^N c {f) iCQ{m 2 x ,rn f )c*lf e[ii,v, ki, k 2 ] f = {q u ,qd,^iJ,xf,xf} (59) 
/ 

where N c (f) is the color factor for the fermions. In the domain < m x < 2m f, which is the relevant 
domain for our study, being the axion very light, the pseudoscalar triangle when both photons are on 
mass-shell is given by the expression 

m ((Am 2 Y l/2 \ m 

Co(m 2 m f ) = ^75- arctan 2 ( I — / - 1 ] I = — t^tV( t ) ( 60 ) 

x 1 7r 2 m x \ \ m 2 , / / 7r 2 m x 
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The coefficient Cy 7 is the factor for the vertex between the axi-Higgs and the fermion current. The 
expressions of these factors are 



V2y u M st u? v 2 



(v 2 v 2 s + v\vl) [4Mf t t; 2 + gBB s (v 2 v 2 s + v\v\ 



V2 Vd 
Vl Vu 



,x><?« 



c x,i = yi c x,id^ 

Vd 

We obtain the following expression for the decay amplitude 

1 



(61) 
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(62) 



where the three terms correspond, respectively, to the point-like WZ term, to the 1-loop contribution and 
to their interference. 

Notice that in the expression of this decay rate both the direct (~ (^ 7 ) 2 ) and the interference (~ g 77 ) 
contributions are suppressed as inverse powers of the Stiickelberg mass, here taken to be equal to 1 TeV. 
We have chosen v as the SM electroweak vev, for vs we have chosen the value of 500 GeV. In order to 
have an acceptable Higgs spectrum, the Yukawa couplings have been set to give the right fermion masses 
of the Standard Model, while for gs and B$ we have chosen gs = 0.1 and B$ = 4. 

We show in Figs. [6] and [7] results obtained from the numerical evaluation of the decay amplitude as 
a function of the mass of the axion m x , which clearly indicates that the decay rates are very small for 
a milli-eV particle, although larger than those of the PQ case 30J. We conclude that a PQ-like axion is 
indeed long-lived also in these models and as such could, in principle, contribute to the relic densities of 
dark matter. For an axion with a mass in the MeV region, instead, the particle is not stable and as such 
would decay rather quickly. The decay, in this case, is fast enough (r < 10~ 3 sec) and does not interfere 
with the nucleosynthesis. 



7 Cold dark matter by misalignment of the axion field 

In the case of a long-lived axion, the generation of relic densities of axion dark matter, in this model, 
involves two (sequential) misalignments, generated, as we have already discussed, the first at the elec- 
troweak scale, and the second at the QCD phase transition. The presence of two misalignments at two 



separate scales, as discussed in 30J], is typical of axions which show interactions both with the weak and 
with the strong sectors, due to the presence of mixed anomalies. This point has been addressed in detail 
within a non-super symmetric model, but in a supersymmetric scenario the physical picture remains the 
same. 
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Figure 6: Decay amplitude (left panel) and mean lifetime (right panel) for \ ~ * 77 as a function of the axi-Higgs 
mass. 



In the PQ-like case, at the first misalignment, taking place at the electroweak scale, the physical axion 
is singled out as a component of Im b, with a mass which is practically zero, due to the small value of the 
curvature induced by the potential generated by electroweak instantons (Fig. [3l top) given in Eq. ([40]) , In 
the case of a very small extra potential (A e jj ~ Knst) it is the second misalignment to be responsible for 
generating an axion mass. At the second misalignment, taking place at the QCD phase transition, the 
mass of this pseudo Nambu-Goldstone mode is redefined from zero to a small but more significant value 
(~ 1CP 3 eV) induced by the QCD instantons ((Fig. [3j bottom). The final value of the mass is determined 
in terms of the hadronic scale Aqcd and of a second intermediate scale, Mg t /v, which replaces f a in all 
of the expressions usually quoted in the literature and held valid for PQ axions, as we are now going to 
clarify. 

• MeV axion 

An MeV axion is allowed only if the extra potential (the misalignment) is assumed to be generated at a 
scale different from the electroweak phase transition, say at an earlier time. This misalignment, in fact, 
should be unrelated to the (quasi-periodic) corrections induced at the electroweak time, as shown in Fig. 
4, being the latter of very small size. However, such an axion would not be long lived. One can easily 
realize that in this scenario, due to the sizeable value of m x , there is an overlap between the period of 
coherent oscillations at the QCD hadron transition and the typical lifetime at which the axion decays. 
This can be trivially checked by comparing the QCD time, defined as the inverse Hubble rate at the 
temperature of confinement (H(Tqcd) ~ 10~ n eV, Tqcd ~ 200 MeV) tQCD ~ 10 -4 sec with the axion 
lifetime in this typical mass range. 

• PQ-like axion 

For a PQ-like axion the effective scale (Mg t /v) is the result of the product of two factors: a first factor 
due to the rotation matrix of the Stiickelberg field Im6 onto x - which is proportional to v/M - times a 
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second factor (1/Mst) which is inherited from the original Im b/M$t FF (WZ) counterterm. Specifically, 
starting from Eq. (|39p . the size of the projection of Imfr into x 1S given by 

-^B s g B (v 2 v 2 s + v\v\) ~ v/Mst (63) 

and hence a typical PQ interaction term involving the Stiickelberg field b becomes 

Imb _x 

f 2 St 



FF' -> £ , FF'. (64) 

M st Mh/v V ; 



The physical state with a b component (i.e. x) acquires an interaction to FF which is suppressed by the 
scale Mgjv. 

Having identified this scale, if we neglect the axion mass generated by V (Eq. I4U|) at the electroweak 
scale, the final mass of the physical axion induced at the QCD scale is controlled by the ratio m x ~ 
AQ CD v/Mg t , where the angle of misalignment is given by 9' = x v /Mg t . 

Coming to the value of the abundances for a PQ-like axion - defined as the number density to entropy 
ratio Y = n x /s - these can be computed in terms of the relevant suppression scale appearing in the xFF 
interaction. We have expanded on the structure of the computation in Appendix [UJ If we indicate with 
0'{Ti) the angle of misalignment at the QCD hadron transition and with Tj the initial temperature at the 
beginning of the oscillations, the expression of the abundances takes the form 

Yx[Ti) = im) v^w^mp' (65) 

which depends linearly on Mst- As we have already mentioned, the computation of the relic densities for 
a non-thermal population follows rather closely the approach outlined in the non-super symmetric case. 



For instance, a rather large value of Mst, of the order of 10 7 GeV 30], determines a sizeable contribution 
of the gauged axion to the relic densities of cold dark matter. These, in turn, follow rather closely the 
behaviour expected in the case of the PQ axion. In practice, to obtain a sizeable non-thermal populations 
of gauged axions, Mst should be such that Mo t /v ~ f a , with f a the usual estimated size of the PQ 
axion decay constant. This allows a sizeable contribution of x to the relic density of cold dark matter, 
with a partial contribution to Q (Q x /i 2 ~ 0.1) in close analogy to what expected in the case of the PQ 
axion. These considerations, which are in close relations with what found in the non-supersymmetric 
construction 30j], in this case will be subject to the constraints coming from the neutralino sector and 
its abundances derived from WMAP. We will come back to this point after presenting the results of our 
simulations in the next sections. 



8 The neutralino sector 

The neutralino sector is constructed from the eigenstates of the space spanned by the neutral fields (i \\y3 , 
i\y, i^B, H\, H^, S, iph)i which involve the three neutral gauginos, the two Higgsinos, the singlino (the 
fermion component of the singlet superfield) and the axino component of the Stiickelberg supermultiplet. 



24 



v s = .5 TeV M st = 1 TeV tanp = 35 



v s = .5 TeV Ms, = 1 TeV tanp = 35 



1e-20 




1e-25 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



0.01 



0.001 



1e-04 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



Figure 7: Decay amplitude and mean lifetime for \ ~~ ^ 11 as a function of the axi-Higgs mass for an axion whose 
mass in the MeV range. 



We denote with M xo the corresponding mass matrix and we list its components in the appendix. The 
neutralino eigenstates of this mass matrix are labelled as Xi (i = 0, . . . , 6) and can be expressed in the 
basis {iXy/ z , iXy, H\, H\, S, ipb} 



an i\w 3 + a i2 iXy + iXs + a i4l Hi + a i5 H% + a i6 S + a i7 ipb- 



(66) 



The neutralino mass eigenstates are ordered in mass and the lightest eigenstate corresponds to i = 0. We 
indicate with O x the rotation matrix that diagonalizes the neutralino mass matrix. In order to perform 
a numerical analysis of the model we need to fix some of the parameters, first of all requiring consistency 
of their choice with the masses of the Standard Model particles. For this purpose, the Higgs vev's v\ 
and V2 have been constrained in order to generate the correct mass values of the W^, which depends on 
v 2 = vf +v%, and of the Z gauge boson. 

The Yukawa couplings have been fixed in order to give the correct masses of the SM fermions. The 
choice of vs and of the parameter A in the trilinear term XSH\ ■ H 2 in the scalar potential has been made 
in order to obtain mass values in the Higgs sector in agreement with the limits from direct searches (with 
A ~ O(l))- For this reason we have selected the value A = 0.5 and the assignment Bh 1 = —l,Bg = 3 for 
the U(1)b charges of the Higgs and the singlet; Bq = 2 for the quark doublet, and Bl = 1 for the lepton 
doublet. The gauge mass terms parameters have been selected according to the relation 



My : M w : Mq = 1:2:6, 



(67) 



coming from the unification condition for the gaugino masses. As a further simplification, the sfermion 
mass parameters Ml, Mq, ttir, mo and mjj have been set to a unique value Mq. We have also chosen 
a common value ao for the trilinear couplings a e ,ad and a u . With these choices, besides tan/3, the only 
other free parameters left are the Stiickelberg mass Mst, the gaugino mass term for A#, denoted by Mb, 
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and the axino mass term, Mb- Our choices are the following 



M Y = 500GeV M YB = 1 TeV M w = 1 TeV M G = 3 TeV 
Ml = Mq = rriR = mo = ra t/ = Mq = 1 TeV 
a e = ad = a u = a = 1 TeV 

a A = -100GeV, (68) 

where with Ml and Mq we have denoted the scalar mass terms for the sleptons and the squarks, assumed 
to be equal for all the 3 generations. We have also chosen 

M B = M b = 1 TeV (69) 

with a coupling constant gs of the anomalous U(l) of 0.4. From previous investigations such values of the 



anomalous coupling are known to be compatible with LEP data at the Z resonance 43|, |44|. In particular, 



the mass of the extra Z' (Mz 1 ), which in our case is of the order of the Stiickelberg mass, Mz> ~ Mst 



431 ] . due to the region of variability of Mst that we investigate in our simulations, obviously satisfies the 



current LHC constraints at 95% CL ( > 1140 GeV) from CMS [45j] and from ATLAS (> 1.83 TeV) [46 1 
on the absence of a resonance in the dilepton channel at 7 TeV for an extra Z prime with Standard-Model 
like couplings. This parameter choice is our benchmark, which is compatible with all the SM requirements 
on the spectrum of the known particles. It involves SUSY breaking scales in a kinematical range which is 
under investigation at the LHC. We also assume a value Vb = 20 GeV for the vev of the saxion field Re6. 

The most significant parameters in the relic density calculation are Mst and the Higgs vev ratio tan j3. 
Concerning the Stiickelberg mass, its value has been chosen to be varied in two different regions, 2 — 10 
TeV and 11 — 25 TeV. In both regions we will consider different values of tan/3. 



9 Neutralino relic densities and cosmological bounds 

As it is well known, the evaluation of the relic densities requires the calculation of a great number of 
thermally averaged cross sections, given the number of particles which are present. Before coming to the 
discussion of the results of this very involved analysis, which is summarized just in some simple plots of 
the relic densities of the lightest neutralino - as a function both of Mst and tan j3, - we present a general 
description of the structure of the interactions in the model. We also list the 2-to-2 processes that have 
been considered in the coupled Boltzmann equations. 

We start from the action involving the physical axion (Hq ) and its interactions with the various sectors. 
These involve, typically, interactions with the Higgs sector via bilinear vertices (proportional to R H o HH ), 
and trilinear ones (proportional to R H o HHH ^j in H, with H denoting generically CP-even and CP-odd 
Higgs eigenstates. Other interactions in the same component of the Lagrangian involve axion-neutralino 
terms (R H ° XiX j) plus axion-charginos (R H ° Xi x i ). Other terms are those involving interactions of the 
axion with the sleptons (R H oh^) and the squarks (R H o^i^)- : vertices involving gauge bosons (for instance 
R H » Ax i x f ; with a photon A and two charginos) and quartic contributions with 2, 3 and 4 axion lines. 



26 




The Lagrangian describing all the tree-level interactions involving the axion is 

C H s_ mt =R$ H t Hi oH*H*lP + R H o 2H » {Hi) 2 ^ + R^^HlHtHZHi + R^^l ( H *) 2 H*H> + 
R H o H o s H 5 H^ + R H ° 2H S 2 {Hl f {Hi) 2 + i^o 53 "o 4 (^) 3 flg + R H o 4 (# 5 ) 4 + 
R H o^Hlxh°j + R Hf ^HlxfxJ + R H ^Hl~l}l 3 + R H ^Hlq\q J + 
R^^lHlA^xHxf + R H ° Zx > x >HlZ»xtwJ + R H ° Z ' x ^HlZ'»xh, X J+ 
R H ° WTx * X °HlW?xyxT + R xh±H ^xhfH^Hl + R^^xh'jHlHl. (70) 



The explicit expressions of these vertices are rather involved and we omit them. Other interactions 
appearing in the interaction Lagrangian involve derivative couplings with the gauge bosons and the 
Higgses and they are given by 

£ H B_ int = R^x^HlW^d^H* + R H ° H » A HlA^Hl + R H o H « z H^^H^ 

R H l H l z ' HlZ'^Hl (71) 

Similar interactions are also typical for H^ the CP-odd Higgs. Some of the vertices are illustrated in 
Fig. [8l Besides the interaction with the axi-Higgs, we have the following vertices involving neutralinos 

Cn-M = R^z^Hx* + R x ^ z ' z'»xh»x* + R x ^ w ^xh,xf + R x ^H^xf+ 
R x^ H k x o x o + sfiff^ffo + R^^xHm + R x ^ I x°xfflfi+ 

R^t^^xfH^H* + R^ H ^xh ^H l + R^ H ^ X l X ]H^. (72) 
Some of these vertices are illustrated in Fig. [9l The full lagrangian has been implemented using the 
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Figure 9: Neutralino interactions 
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Table 1: Tree level neutralino annihilation processes in the 3 kinematic channels 



FevnRulesj47f| package. The same package allows to generate the CalcHEP 48] model files which are 



needed by micrOMEGAsH 

for the calculation of the scattering cross section that are required in the 

relic density calculation. 

With our choice for the parameters the lightest neutralino is the Lightest Supersymmetric Particle 
(LSP) and so it is the dark matter component in our simulations. The value of the neutralino mass in this 
case turns out to be around 23 GeV with a rather mild dependence on tan f3. For tan /3 varying between 
5 and 25 the neutralino mass varies from 22.4 to 23.8 GeV. 

We show in Tab. [1] a list of the most relevant 2-to-2 processes which are generated in the s,t and 
u channels having neutralinos in the initial state (in), while the possible final states are shown on the 
right-hand side of the same table (out). 

In Fig. [10] we show the results obtained for the lightest neutralino relic density with M$t in the range 
5 — 8 TeV, vs = 600 GeV and a varying tan f3. The values of vs, tan f3 and Mst f° r which we plot the result 
coming from the relic density calculation are those that give also acceptable mass values for the whole 



spectrum, in particular for the neutral Higgs (~ 124 — 126 GeV) 50. l51|. The horizontal bar represents 
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Figure 10: Relic density of the lightest neutralino as a function of the Stiickelberg mass in the region 5 
for different values of the Higgs vev ratio tan j3 



8.5 TeV 



the experimental value for the physical dark matter density measured by WMAP, Clh 2 = . 1 123 ± .0035 52]. 
In Fig. [11] we show the analogous results obtained in the range 11 — 24 TeV with vs = 1-2 TeV and varying 
tan (3. Once again these values are such that we obtain acceptable values for the masses of all the particles 
in the model. One can immediately notice that for a fixed value of tan/3 as we increase Mst, the relic 
densities grow and tend to violate the WMAP bound. This trend has been found over a sizable range of 
variability of tan /? and is a central feature of the model. It is then obvious, from the same figures, that 
it is possible to raise the Stuckelberg mass and stay below the bound if, at the same time, we increase 
tan (3. 



10 Summary: windows on the axion mass 

At this point, before coming to our conclusions, we can try to gather all the information that we have 
obtained so far in the previous sections, summarizing the basic properties of axions in these types of 
models. 

• The milli-eV (PQ-like) axion 

One possibility that we have explored in this work is that V, the extra potential which is periodic in the 
axion field, may be generated around the TeV scale or at the electroweak phase transition. The actual 
strength of the potential, remains, in our construction, undetermined and the physical features of the axion 
(primarily its mass), depend upon this parameter. We have tried to describe the various possibilities, in 
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Figure 11: Relic density of the lightest neutralino as a function of the Stuckelberg mass in the region 2 — 24 TeV 
for different values of the Higgs vev ratio tan /3 



this respect, and the essential features for each choice for the value of the mass. In particular, if the extra 
potential is generated by non-perturbative effects at the electroweak phase transition, then the mass of 
the axion is tiny and the true mechanism of misalignment which determines its mass takes place at a 
second stage, at the QCD phase transition. In this case the physical axion of the model would be no much 
different from an ordinary PQ axion and would be rather long-lived. At the same time, its abundances 
are fixed by the possible value of the scale M$ t /v, which should be rather large (~ 10 10 — 10 12 GeV), 
of the same order of f a in typical axion models, to be a significant component of cold dark matter. In 
the region that we have analyzed numerically, with M$t around the 2-20 TeV's, the contribution to dark 
matter from misalignment of the axion field, in this case, should be small. 

A second important constraint on this particle, in this mass range, comes from direct axion searches, 
which also requires the interaction of the axions with the gauge fields (in particular the photon) to be 
suppressed by a large f a . For this reason, with M$t in the TeV region, these simulations indicate that an 
axion of this mass, in fact, can be excluded by typical searches with detectors of Sikivie type. The reason 
is rather obvious, since axions in the milli-eV mass range could be copiously produced at the center of the 
sun and probably should have been seen by now in ground based detectors (helioscopes), such as CAST 
5a] . We recall that one of the goals of searches with helioscopes is to set a lower bound on the suppression 
scale f a of the axion-photon vertex, which is currently experimentally constrained, as we have already 
mentioned, to be rather large. 



The meV axion 
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A second possibility that we have investigated is that the extra potential appearing in the CP-odd sector is 
unrelated to instanton corrections in the electroweak vacuum. In this case the mass of the axion remains 
a free parameter. The range that we have explored in this second case involves an axion mass in the 
MeV region, discussing the several constraints that emerge from the model. In this case the axion is, in 
general, not long-lived and as such is not a component of dark matter. On the other hand, the constraints 
from CAST can be avoided, since the particle would not be produced by ordinary thermal mechanisms 
at the center of the sun, where the temperature is about 1.7 keV, being its mass above the keV range. 
Obviously, in this case other constraints emerge from nucleosynthesis requirements, since a particle in this 
mass range has to decay fast enough in order not to generate a late entropy release at nucleosynthesis 
time. We have seen that an axion in the MeV range is consistent with these two requirements. An axion 
of this type could be searched for at colliders, and in this respect the analysis of its possible detection 
at the LHC would follow quite closely the patterns described by two of us in 4l|. As in this previous 
(non-supersymmetric) study, where the axion is Higgs-like (of a mass in the GeV region) typical channels 
where to look for this particle would be a) the associated production of an axion and a direct photon; 
b) the multi-axion production channel, and c) the associated production of one axion and other Higgses 
of the CP-even sector. The modifications, compared to that previous study, would now involve 1) the 
lower value of the mass of the axion (MeV rather than GeV); 2) the presence of extra supersymmetric 
interactions. 



10.1 Comments 

One important comment concerns the connection between these class of models and their completion 
theories such as string theory, which lay at their foundation. In our study we have selected a scenario 
characterized by low energy supersymmetry, with a phenomenological analysis that is essentially connected 
with the TeV scale and above. This is the scale which is likely to be scanned in the near future by several 
experiments, including the LHC, and for this reason we have directed out numerical studies in this 
direction. There is, however, a second it involves a value of M$t which is very large and close to the 
Planck scale. In this second case the model predicts, obviously, a decoupling of the anomalous symmetry, 
leaving at low energy a scenario which is essentially the same of the MSSM, since the extra Z prime, 
which is part of the spectrum, is extremely heavy. This would obviously imply a decoupling both of 
the anomalous gauge boson and of the anomalous trilinear interactions which are associated with it. A 
physical axion could, however, survive this limit, if the scale of the extra potential is also of the order of 
Mst, but its interaction with ordinary matter would be extremely suppressed by the same scale. 



11 Conclusions 

The investigation of the phenomenological role played by models containing anomalous gauge interactions 
from abelian extensions of the Standard Model, we believe that will receive further attention in the 
future. These studies can be motivated within several scenarios, including string and supergravity theories, 
in which gauged axionic symmetries are introduced for anomaly cancellation. In turn, these modified 
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mechanisms of cancellation of the anomalies, which involve an anomalous fermion spectrum and an axion, 
are essentially connected with the UV completion of these field theories, which in a string framework is 
realized by the Green-Schwarz mechanism. 

The model that we have investigated (the USSM-A) summarizes the most salient physical features of 
these types of constructions, where a Stiickelberg supermultiplet is associated to an anomalous abelian 
structure in order to restore the gauge invariance of the anomalous effective action. In this work we 
have tried to characterize in detail some among the main phenomenological implications of these models, 
which are particularly interesting for cosmology. The physical axion of this construction, or gauged axion, 
emerges as a component of the Stiickelberg field Im b. We have pointed out that the mechanism of sequen- 
tial misalignment, formerly discussed in the non-supersymmetric case [30] , finds a natural application also 
in the presence of super symmetry, with minor modifications. 

One relevant feature of these models, already noticed in [3(3], is that their axions do not contribute 
to the isocurvature perturbations of the early universe, being gauge degrees of freedom at the scale of 
inflation. 

We have followed a specific pattern in order to come out with specific results in these types of models, 
using for this purpose a particular superpotential (the USSM superpotential), whose essential features, 
however, may well be generic. 

We have presented an accurate study of the neutralino relic densities, showing that the Stiickelberg 
mass value is constrained by the requirement of a consistent mass spectrum, with values for the lightest 
CP-even Higgs larger than the current LHC limits ( > 120 GeV) 

and by the experimental value for the dark matter abundance from WMAP 52J. Thus, in these 
models, the allowed value of the Stiickelberg scale is positively correlated with the value of tan (3. As 
it grows, tan (3 has also to grow (for a fixed value of the vev of the singlet vs) in order to preserve the 
WMAP bound. In particular M$t and vs are positively correlated. This correlation is necessary in order 
to obtain values of the neutralino mass which allow to satisfy the same bounds, which in our case is 
around 20 GeV. 

We have seen that with a Stiickelberg mass in the TeV range the non-thermal population of axions does 
not contribute significantly to the dark matter densities if these axions are PQ-like. These types of con- 
straints, obviously, are typical of super symmetric constructions and are avoided in a non-supersymmetric 
context. In this second case, as discussed in 3CJ], a Stiickelberg scale around 10 7 GeV is sufficient to revert 



this trend. 

We have also pointed out that gauged axions in the milli-eV mass range are probably difficult to 
reconcile with current bounds from direct searches, while the case for detecting MeV or heavier axions, in 
these types of models, remains a wide open possibility. In this second case, cascade decays of these light 
particles and their associated production with photons should be seen as their possible event signatures 
at the LHC. 
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A General Features of the Model 



In this appendix we summarize some of the basic features of the USSM-A. The gauge structure of the 
model is of the form SU(3) C x SU (2) x U (l)y x U(1)b, where B is the anomalous gauge boson, and with a 
matter content given by the usual generations of the Standard Model (SM). In all the Lagrangians below 
we implicitly sum over the three fermion generations. A list of the fundamental superfields and charge 
assignments is summarized in Tab. EJ The Lagrangian can be expressed as 



&USSM-A = C-USSM + C-KM + £-FI + £>a 



(73) 



where the Lagrangian of the USSM (Cussm) has been modified by the addition of C ax ion to compensate 
for the anomalous variation of the corresponding effective action due to the anomalous charge assignments. 
The former is given by 



£-USSM = £>lep + £quark + ^Higgs + C-gauge + C-SMT + &GMT 



(74) 



with contributions from the leptons, quarks and Higgs plus gauge kinetic terms. The matter contributions 
from leptons and quarks 



c, 



quark 



d A e 



d 4 



li e 2g2W+g Y Y+g B Bl + £f ^W+gyY+ggB £ 



(75) 



q] e 2g s G+2g 2 W+g Y Y+g B B Q + (jl e 2g s G+g Y Y+g B Bfj R + 2g s G+g Y Y+9BB f) R (jty 



are accompanied by a sector which involves two Higgs SU(2) doublet superfields, Hi and H 2 , and one 
singlet S 



tgg» 



d 1 



jj^W+gyY+gsB^ + ^ e 2g 2 W+g Y Y+g B B ^ + & e g B Bg + W< 5 2 (0) + W5 2 {9)] (77) 



with the superpotential chosen of the form 

W = XSH 1 • H 2 + y^ • LR + y d H x ■ QD R + y u H 2 ■ QU R . 



(78) 



This superpotential, as shown in 34], [35j, allows a physical axion in the spectrum. The gauge content 
plus the soft breaking terms in the form of scalar mass terms (SMT) are identical to those of the USSM 

Cgauge = \ J d 4 9 [Q a G a + W a W a + W Ya W% + W Ba w£] 5 2 {9) + h.c. 

Csmt = - J d 4 9 5 A {9, 9) [M 2 L DL + m 2 R &R + MqQ^Q + mf;U R U R + m 2 D D R D R 

+ rr^HlHi + m\R\H 2 + m 2 s &S + (a A S#i • H 2 + h.c.) + (a e H 1 • LR + h.c.) 

+ {a d H 1 • QDr + h.c.) + (a u H 2 ■ QU R + h.c.)}. (79) 
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Table 2: Charge assignment of the model 



As usual, Ml, MQ,mR,mu R ,TnE) R ,rni,rn2,rns are the mass parameters of the explicit supersymmetry 
breaking, while a e ,a\, a u >&d are dimensionful coefficients. The soft breaking due to gaugino mass terms 
(GMT) now include a mixing mass parameter Myb 



£gmt = I d ft 



\ {M G g a g a + M w W a W a + M Y W Ya Wl + M B W Ba W B 
+M Y bW Yq W b ) + h.c] 5 A {9, 9). (80) 



The superfield b describes the Stiickelberg multiplet, 

b = b + V29if) b - iOa^ed^b + -±=0005*8^ ~ \o999Ub - 99F b , (81) 

v2 4 

and contains the Stiickelberg axion (a complex b field) and its super symmetric partner, referred to as the 
axino (ipb), which combines with the neutral gauginos and higgsinos to generate the neutralinos of the 
model. Details on the notation for the superfields components can be found in Tab. We just recall that 
we denote with X B and Ay the two gauginos of the two vector superfields (B,Y) corresponding to the 
anomalous U(1)b and to the hypercharge vector multiplet. The singlet superfield S has as components 
the scalar "singlet" S and its super symmetric partner, the singlino, denoted as S. 

The interactions and dynamics of the axion superfield are defined in C ax ion, the Lagrangian that 
contains both the kinetic (Stiickelberg) term, responsible for the mass of the anomalous gauge boson 
(which reaches the electroweak symmetry breaking scale already in a massive state), the kinetic term of 
the saxion and of the axino, and the Wess-Zumino terms, which are needed for anomaly cancellation. We 
recall that Stiickelberg fields appear both in anomalous and in non-anomalous contexts. The second one 
has been analyzed recently in [541 ]. 

The extra contributions to C a xim, called C ax ion,i are given by 

£axion,i =M st RebD B + M st (iipb^B + h.c.)+ 
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Table 3: Superfields and their components. 

F h X G X% + l -^G a u ^ h ^a u \ G + — ^Im6 G atxv G a av + 

M st \l6 b G G 8V2 G G ^ 16v/2 G 64^/2 

-^Im& XqOuD^Xq + -^Re6 D G D G ^Re6 G a ^G" - ^Reb X a G aaD^X a G 

8^/2 G M G 16^2 16^2 M 8\/2 G ^ G 
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^=Im6 X^D^ + ^Refo D^ w - ^Reft ^ v\ 

+ 7!^ A ^ b " ^^"^ + ^ 

^Re5 DyDy ^Re6 F r ^ - ^Reb Xya^Xy + h.c) + 

2\/2 2\/2 v2 ' 

^ (\F h X B X B + -^=D_BA_BVb - ^=<,Vb^^A B + ^=Im& F B - -±=Imb X B a^X B + 
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8^2' 

We finally recall that the three scalar sectors of the model are 



characterized in terms of 



(82) 
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• A Charged Higgs sector 

This sector involves the states (ReH^ , ReHf). The mass matrix has one zero eigenvalue correspond- 
ing to a charged Goldstone boson and a mass eigenvalue corresponding to the charged Higgs mass 

2 fvi V 2 \ (1 2 1 2 , V S \ (oo\ 

m H±=[- + ~) {- A 9 v lV2 - -A v lV2 + a x -j=J ; (83) 

where g 2 = g%+9y- 

• A CP-even sector 

This sector is diagonalized starting from the basis (Rei?|, ReiJf , ReS*, Re6) The four physical states 
obtained in this sector are denoted, as Hq,Hq,Hq and Hq. Together with the charged physical 
state extracted before, H^, they describe the 6 degrees of freedom of the CP-even sector. 

• A CP-odd sector 

This sector is diagonalized starting from the basis (Imii^ , ImH 2 j hnS 1 , Im6) . We obtain two physical 
states, Hq and Hq, and two Goldstone modes that provide the longitudinal degrees of freedom for 
the neutral gauge bosons, Z and Z'. 



B Neutralino mass matrix 



Now we turn to the neutralino sector; the mass matrix in the basis (i\ W3 ,i\y ,iX B , H\, H 2 , S,ipb) takes 
the form 



with 



/Ml 1 
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M 22 
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M 2 f 
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M x 3 3 
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M 35 
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(84) 
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M x» 2 
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9^1 M 15 _ 52^2 17 _ n 
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1 M M 24 9YV1 25 ffT«2 
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The rotation matrix for this sector is implicitly defined as O x and 



f i^w 3 ' 




(x°o) 


iXy 




X? 


iX B 






H\ 




xl 


Hi 




xl 


S 




xl 






\xV 



(86) 



B.l Chargino sector 

We recall here the structure of the chargino sector and the diagonalization procedure. We define 

X w + = — j={X Wl iX W2 ) = ~^J^-^ Wl 

and in the basis (^X w +, H\ , A„,- , Hf^j we obtain the mass matrix 

/ M w g 2 v 1 \ 

g 2 v 2 Xv,s 

M w g 2 v 2 

\g 2 Vi Xv s J 

from the diagonalization we get the squared eigenvalues 



(87) 



M| ± 



(88) 



m Y ± = O 



M w + ^4 + gW T \J (M w + X*v 2 s + g*v 2 ) 2 _ 4 ( Aus M w - si^y 



If we define 



Hi 



(89) 
(90) 



and define the mass eigenstates as 

X + = Vi(j + X ~ = Ui/j- (91) 
where U and V are two unitary matrices that perform the diagonalization of this sector. If we define 



X 



M w g 2 v 2 
a 2 vi Xv s , 



then these unitary matrices are defined is such a way that 

VX^XV- 1 = U*XX^U T = M x±4iag - 

where M x ± >diag is given by 



l X ± ,diag 



m~± 
m~± 

A2 



(92) 



(93) 



(94) 
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C Appendix. Relic densities at the second misalignment 



In this appendix we fill in the gaps in the derivation of the expression of the abundances generated by 
the mechanism of vacuum misalignment. We start from the Lagrangian 

S = f d A x^ Qx 2 - \rn 2 x Y x ±\ , (95) 

where T x is the decay rate of the axion and we have expanded the potential around its minimum up to 
quadratic terms. The same action is derived from the quadratic approximation to the general expression 

S = J d A xR 3 (t) QoJ (d a 6) 2 - p 4 (1 - cos d) - V^j (96) 

which, in our case, is constructed from the expression of V given in Eq. (133p . with p ~ v, the electroweak 
scale. We also set other contributions to the vacuum potential to vanish (Vo = 0). In a Friedmann- 
Robertson- Walker spacetime metric with a scaling factor R(t), this action gives the equation of motion 



d_ 

li 



{R 3 (t)(x + r, 



+ R 3 m 2 (T) = 0. (97) 



We will neglect the decay rate of the axion in this case and set T x ~ 0. At this point, since the potential 
V is of non- perturb at ive origin, we can assume that it vanishes far above the electroweak scale (or 
temperature T ew ). For this reason m x = = for T 3> T ew , which is essentially equivalent to assume 
that the Stiickelberg axion is not subject to any mixing far above the weak scale. The general equation 
of motion derived from Eq. (197p . introducing a temperature dependent mass, can be written as 

X + 3Hx + m x (T) X = 0, (98) 

which clearly allows as a solution a constant value of the misalignment angle = 6%. The T-dependence of 
the mass term should be generated, for consistency, from a generalization to finite temperature of V'. In 
practice this is not necessary in our case, being the role of the first misalignment negligible in determining 
the final mass of the axion. 



The axion energy density is given by 



which after a harmonic averaging gives 



1 .o 1 



P=^X 2 + ^ 2 X X\ (99) 



(p) = m 2 x ( X 2 ). (100) 

Notice that after differentiating Eq. (|99p and using the equation of motion in (|98p . followed by the 
averaging Eq. (jlOOp one obtains the relation 

(p) = (p) (-3H + -) , (101) 
V ml 
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where the time dependence of the mass is through its temperature T(i), while H(t) = R(t)/R(t) is the 
Hubble parameter. By inspection one easily finds that the solution of this equation is of the form 

showing a dilution of the energy density with an increasing space volume, valid even for a T-dependent 
mass. At this point, the universe must be (at least) as old as the required period of oscillation in order 
for the axion field to start oscillating and to appear as dark matter, otherwise 9 is misaligned but frozen; 
this is the content of the condition 

m x (Ti) = 3H(Ti), (103) 

which allows to identify the initial temperature of the coherent oscillation of the axion field x-> by 
equating m x (T) to the Hubble rate, taken as a function of temperature. 

To quantify the relic densities at the current temperature To (To = T(to), at current time to) we define 
preliminarily the two standard effective couplings 

3 - ,m \ 3 



B F 
4 



fc,T = Eft('£Y+^£ft(!) . (104) 



B 



functions of the massless relativistic degrees of freedom of the primordial state, with T 3> T ew . The 
counting of the degrees of freedom is: 2 for a Majorana fermion and for a massless gauge boson, 3 for 
a massive gauge boson and 1 for a real scalar. In the radiation era, the thermodynamics of all the 
components of the primordial state is entirely determined by the temperature T, being the system at 
equilibrium. We exclude for simplicity all sorts of possible source of entropy due to any inhomogeneity 



(see, for instance, 55]). Pressure and entropy are then just given as a function of the temperature 



7T 2 

P = 3p = — g, tT T 
2vr 2 

a = -tt9*,s,tT\ (105) 
45 

while the Friedmann equation allows to relate the Hubble parameter and the energy density 



H = ^-nG NP , (106) 

with Gn = 1/Mp being the Newton constant and Mp the Planck mass. The number density of axions 
n x decreases as 1/R 3 with the expansion, as does the entropy density s = S/R 3 , where S indicates 
the comoving entropy density - which remains constant in time (5 = 0)- leaving the ratio Y a = n x /s 
conserved. We define, as usual, the abundance variable of x 



(107) 

Ti 
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at the temperature of oscillation Tj, and observe that at the beginning of the oscillations the total energy 
density is just the potential one 

p x = n^m^Ti) = l/2mj(ri)x?. (108) 
We then obtain for the initial abundance at T = Ti 

2 s 47r z g^ >s ,TTf 

where we have inserted at the last stage the expression of the entropy of the system at the temperature 
Ti given by Eq. (|105j) . At this point, plugging the expression of p given in Eq. (|105p into the expression of 
the Hubble rate as a function of density given n Eq. (|106p . the condition for oscillation Eq. f|103j) allows 
to express the axion mass at T = Ti in terms of the effective massless degrees of freedom evaluated at the 
same temperature, that is 

This gives for Eq. (|109j) the expression 

YJTi) = — x 1 , (111) 

where we have expressed x m terms of the angle of misalignment 6i at the temperature when oscillations 
start. We assume that Oi = (6) is the zero mode of the initial misalignment angle after an averaging. As 
we have already mentioned, Ti should be determined consistently by Eq. (]103p . However, the presence of 
two significant and unknown variables in the expression of m x , which are the coupling of the anomalous 
U(l), gs, and the Stiickelberg mass M, forces us to consider the analysis of the T-dependence of x 
phenomenologically less relevant. It is more so if the Stiickelberg mass is somehow close to the TeV 
region, in which case the zero temperature axion mass m x acquires corrections proportional to the bare 
coupling (ra x ~ \v(\ + 0(<?b)). 

For this reason, assuming that the oscillation temperature Tj is close to the electroweak temperature 
T ew , Eq. (jllOp provides an upper bound for the mass of the axion at which the oscillations occur, assuming 
that they start around the electroweak phase transition. Stated differently, mass values of x such that 
m(Tj) <C 3H(Ti) correspond to frozen degrees of freedom of the axion at the electroweak scale. This is 
clearly an approximation, but it allows to define the oscillation mass in terms of the Hubble parameter 
for each given temperature. 

We recall that the relic density due to misalignment can be extracted from the relations 

n mis _ P%1_ = i.n x0 m x ) = / n^o\ m x s 

X Pc Pc \ S J Pc 

where we have denoted with n x o the current number density of axions and with p™§ s their current energy 
density due to vacuum misalignment. This expression can be rewritten as 



X s 



(113) 
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using the conservation of the abundance Y a Q = Y a (Ti). Notice that in Eq. (I113P we have neglected a 
possible dilution factor 7 = s osc / sq which may be present due to entropy release. We have introduced the 
variable 

- - 3H ° (114) 



r 8ttG n 
which is the critical density and 

2?r 2 

S ° = -J^9*s> T o T o ( 115 ) 
which is the current entropy density. To fix g*s,T we J us t recall that at the current temperature T the 
relativistic species contributing to the entropy density sq are the photons and three families of neutrinos 
with 

7 o nf T A 3 



5*5,T =2 + -X3x2l^l (116) 

where, from entropy considerations, T v /Tq = (4/11) 1 / 3 . 

To proceed with the computation of the massless degrees of freedom above the electroweak phase 
transition we just recall the structure of the model. We have 13 gauge bosons corresponding to the gauge 
group SU(3) x SU(2) x Uy (1) x [7^(1), 2 Higgs doublets, 3 generations of leptons and 3 families of quarks. 
Above the energy of the electroweak transition we have only massless fields with the exception of the 
Ub(1) gauge boson, since this symmetry takes the Stiickelberg form above the electroweak scale, giving 
g^^T = 110.75. Below the same scale this number is similarly computed with g* y = 91.25. Other useful 
parameters are the critical density and the current entropy 

p c = 5.2 • 10~ 6 GeV/cm 3 s = 2970 cm~ 3 , (117) 

with 9i ~ 1. It is clear, by inserting these numbers into Eq. (|112|) that Q r y 



2y/$mg %Ti TiMp Pc 

is negligible unless a x ~ Mo t /v is of the same order of f a ~ 10 12 GeV, the standard PQ constant. This 
choice would correspond to U x ~ 0.1, but the value of M$t should be of O(10 ) GeV . 
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